Decaying Dark Matter as a Probe of Unification and TeV Spectroscopy 
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In supersymmetric unified theories the dark matter particle can decay, just like the proton, through 
grand unified interactions with a lifetime of order of ~ 10 26 sec. Its decay products can be detected 
by several experiments - including Fermi, HESS, PAMELA, ATIC, and IceCube - opening our first 
direct window to physics at the TeV scale and simultaneously at the unification scale ~ 10 GeV. 
We consider possibilities for explaining the electron/positron spectra observed by HESS, PAMELA, 
and ATIC, and the resulting predictions for the gamma-ray, electron/positron, and neutrino spectra 
as will be measured, for example, by Fermi and IceCube. The discovery of an isotropic, hard 
gamma ray spectral feature at Fermi would be strong evidence for dark matter and would disfavor 
astrophysical sources such as pulsars. Substructure in the cosmic ray spectra probes the spectroscopy 
of new TeV-mass particles. For example, a preponderance of electrons in the final state can result 
from the lightness of selectrons relative to squarks. Decaying dark matter acts as a sparticle injector 
with an energy reach potentially higher than the LHC. The resulting cosmic ray flux depends only 
on the values of the weak and unification scales. 



INTRODUCTION 

In grand unified theories the proton can decay be- 
cause the global baryon-number symmetry of the low en- 
ergy Standard Model is broken by physics at the grand 
unification (GUT) scale. Indeed, only local symmetries 
can guarantee that a particle remains exactly stable, 
whereas global symmetries are generally broken in fun- 
damental theories. Just as the proton is long-lived but 
may ultimately decay, other particles, for example the 
dark matter, may decay with long lifetimes. It would 
seem miraculous that the dark matter particle's lifetime 
is in a range which is long enough to be a good dark mat- 
ter candidate yet short enough to be observable today. 
Nevertheless, this may well be what happens in Grand 
Unified Theories. For example if a TeV mass dark mat- 
ter (DM) particle decays via GUT-suppressed dimension 
6 operators, its lifetime would be 

' - = 3 x 10- s (™X f „ St V 

ml M Wdm/ V2xl0 16 GeV; 

(1) 

where Mqut ~ 2 x 10 16 GeV is the supersymmetric 
(SUSY) unification scale. 

This lifetime is being probed by several cur- 
rent experiments, as shown in Table I. This can 
be understood, at least for the satellite and balloon 
experiments, because these all generally have simi- 



lar acceptances of <~ (1 m 2 )(l yr)(l sr) w 3 x 
10 11 cm 2 s sr. For comparison, the number of inci- 
dent particles from decaying dark matter with a lifetime 
of 1(F s is ~ J 10 kpc £(0.3 jjJj^XlO-* s" 1 ) » 
10 9 cm 2 s 1 sr x , where these could be pho- 
tons, positrons or antiprotons for example, depend- 
ing on what is produced in the decay. This im- 
plies such experiments observe <~ (3 x 10 11 cm 2 s sr) x 
(l0 -9 cm~ 2 s~ 4 sr -1 ) w 300 events. This coincidence 
may allow these experiments to probe physics at the 
GUT scale, much as the decay of the proton and a study 
of its branching ratios would. In fact, HESS, PAMELA, 
and ATIC may have preliminary indications of dark mat- 
ter from the cosmic ray electron/positron spectrum. The 
Fermi satellite will test this by measuring both the elec- 
tron/positron and gamma-ray spectra with significantly 
improved precision. 

Previously [1], we discussed the general framework 
of dark matter decays induced by GUT scale physics and 
its signals at PAMELA, ATIC and Fermi. In this paper, 
we consider the implications of the recently published 
HESS data [2] on this framework and discuss models that 
fit the overall electron/positron spectrum observed by 
PAMELA, ATIC and HESS. We focus on the correlated 
photon and neutrino signals that could be observed at 
Fermi and IceCube respectively. 
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TABLE I: A lower limit on the lifetime of a dark matter 
particle with mass in the range 100 GeV < tuum < 10 TeV, 
decaying to the products listed in the left column. The ex- 
periment and the observed particle being used to set the limit 
are listed in the right column. All the limits are only approx- 
imate. Generally conservative assumptions were made and 
there are many details and caveats as described in [1]. 



THEORETICAL SETUP 

To study the observational consequences of decay- 
ing dark matter in SUSY GUTs one may follow an ef- 
fective field theory approach and consider an extended 
MSSM with higher dimensional operators parametriz- 
ing GUT effects and leading to dark matter decay. A 
detailed analysis of possible higher dimensional opera- 
tors and the ways to generate them from concrete mi- 
croscopic SUSY GUTs was presented in [1]. Here, for 
defmitencss, we will work in the context of the SO (10) 
models described in [1]. As an example, in addition to 
the standard MSSM interactions, we introduce an addi- 
tional vectorlike (16 m , 16 m ) multiplet at the TeV scale 
and IOqut multiplet at the GUT scale. The relevant 
superpotential interactions involving these fields are 

W = A16 m 16 / 10 GU T + ml6 m f6 m + M GUT 10 GU TlOGUT 

(2) 

We will assume that the singlet S rn is the lightest compo- 
nent of the 16 m and will therefore be dark matter. After 
GUT scale matter and gauge fields are integrated out 
one obtains the dimension 5 operator 16 m 16 m 16/16/ in 
the superpotential and dimension 5 and 6 Kahler terms 
16 m 16 m 10^, WlnWmW^Wf involving m-fields. Of all 
these, the only operator that involves two singlet S m 
components of 16 m is the dimension 6 Kahler term yield- 
ing S^Sm 16^16/ (assuming right-handed neutrinos are 
heavy). Consequently, in this model a thermal relic 
abundance of singlet fields is produced through dimen- 
sion 5 decays of the charged components of 16 m close to 
the BBN epoch. These decays are interesting in their 
own right, as they may explain the observed Lithium 
abundances [3]. On the other hand, dimension 6 decays 
between different components of the singlet supcrmulti- 



plct may lead to observable astrophysical signals that we 
discuss in the rest of the paper. 

Note that these decays may go through operators 
generated by integrating out the heavy U(1)b~l gauge 
boson, or by integrating out heavy IOqut fields. In the 
former case decays are flavor universal, while the latter 
gcncrically lead to flavor non-universal decays. In the 
case of flavor non-universal decays, since the decay rate 
scales as the fourth power of the coupling, it is easy to 
have decays to one flavor dominate over the rest. De- 
pending on the relative strength of gauge and superpo- 
tential couplings and the masses of the heavy fields, both 
possibilities can be realized. 

One may worry that this picture could be spoiled by 
lower dimension operators, such as Kahler kinetic mix- 
ings lOgTjrplO/j and 161^16/. However, these are forbid- 
den by R-parity (under which 16 m is even, and IOqut 
is odd), and m-parity under which both 16 m and IOqut 
are odd. 

For simplicity, in this paper we will focus on the case 
in which the scalar s receives a TeV scale vev. In this 
case dimension 6 operators lead to two body decays of 
the singlet fields to the MSSM fields. We are thus lead 
to two interesting observations: 

• In this case dark matter decay products necessar- 
ily contain MSSM superpartners, because direct 
decays of a scalar into two light fermions are sup- 
pressed by helicity. 

• The production of superpartners, combined with 
the generic expectation that sleptons arc lighter 
than squarks leads to decays dominantly into lep- 
tonic channels due to kinematics. 

These lead to a possible connection between the branch- 
ing fraction of dark matter and the spectrum of its de- 
cay products on the one hand, and the supersymmetric 
spectrum and the decay cascades of superpartners on the 
other. 



ASTROPHYSICAL SIGNALS 
Electrons and Positrons 

GUT induced dark matter decays lead to several 
generic expectations for electron/positron spectra. As 
discussed in the previous section, the dark matter is a 
combination of the scalar (s) and fermion (s) compo- 
nents of the S m superfield. The two body decay of dark 
matter will involve sleptons (I, the superpartner of a lep- 
ton) in the final state. The slepton then further decays 



to its partner lepton and the lightest supersymmetric 
particle (LSP), I — > I + LSP, leading to an injection 
spectrum of the lepton that is flat between a lower and 
upper edge [1] . When dark matter is much heavier than 
the superpartners the lower edge is typically a few GeV 
and the upper edge is roughly at 
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The sensitivity of the upper edge in the injection spec- 
trum to the masses of superpartners may lead to inter- 
esting cross checks at the LHC. 

Current data for the positron fraction from 
PAMELA [4] and the e~ + e+ spect rum from the 
ATIC [5] and HESS [2] experiments include interesting 
hints of a possible excess above background. Due to sys- 
tematic uncertainties it is premature to interpret spec- 
tral shapes observed by ATIC and HESS, particularly 
above 100 GeV. However, this flux will be measured in 
the very near future by both Fermi and PAMELA with 
improved statistics and systematics. The HESS experi- 
ment will also extend its measurement to lower energies, 
overlapping with the excess observed by ATIC. It is thus 
interesting to consider several scenarios for the outcome 
of upcoming experiments and their implications in the 
context of our framework. 

For concreteness we will focus on two scenarios of 
possible electron+positron spectra which are roughly 
consistent with current data, given the systematic un- 
certainties, but will be probed further soon: 

Scenario 1 - A smooth HESS signal - The low 
energy data from ATIC (below ~ 75 GeV) is used to esti- 
mate the astrophysical background. We then find decay 
channels to explain the higher energy HESS excess. Two 
simple examples that realize this are shown in figure 1. 

One possibility is that a heavy scalar (of order 6 
TeV or more) is decaying to a pair of smuons which each 
further decay to a muon and a ncutralino. The muons 
decay further to produce electrons and positrons. For 
example, a 8 TeV scalar decaying to 200 GeV smuons 
which then decay to 100 GeV neutralinos is shown in 
figure 1. The spectral shape (see equation (3)) produced 
by such a cascade fits the spectrum observed by HESS 
quite well. 

A second possibility shown in figure 1 is that a heavy 
6 TeV fermion decays to tau+stau, producing a similar 
spectral shape. A decay into the third family may be 
motivated by minimal flavor violation. In this case the 
observed flux is dominated by the electrons and positrons 
produced in the subsequent decay of the hard (~ 3 TeV) 




FIG. 1: The electron+positron spectrum produced by the 
decay of dark matter (s or s) to smuon pairs (solid black) 
and tau-stau pairs (dot-dashed red) as discussed in the text. 
The HESS and ATIC data is shown by red squares and blus 
circles respectively. The systematic error of HESS is shown 
as a grey band and its superimposed energy uncertainty is 
shown as a lighter grey band. The background and smuon 
signal components of the flux are shown by dotted lines. 



tau[41]. The stau will produce an additional soft com- 
ponent to the spectrum which is subdominant. 

In this scenario, we take the lifetimes of the dark 
matter to be 10 26 seconds for the smuon case and 6 x 10 25 
seconds for the tau case. Though both of these possibil- 
ities produce a similar electron plus positron spectrum 
they may be distinguishable by gamma ray and neu- 
trino observations, as discussed later in this paper. It 
is interesting to note that both of these possibilities re- 
quire flavor non-universal decays of dark matter which 
may naturally be produced in the model discussed above. 
Furthermore, since the decay rate scales as the fourth 
power of the coupling, these flavor non-universal decays 
can be easily dominated by decays to one flavor over the 
others. 

Scenario 2 - Multiple features - We now focus 
on some interesting spectra that consist of multiple fea- 
tures and that occur in simple scenarios of GUT induced 
dark matter decays. We will not necesarily assume that 
either the HESS or ATIC spectra are correct, but rather 
pick three examples to demonstrate some of the generic 
possibilities. 

In the first example, dark matter decays into lep- 
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FIG. 2: The "multi feature" electron+positron spectra pro- 
duced in simple dark matter decays scenarios. The short- 
dashed purple line is an example of multiple features arising 
from decays into leptons and sleptons. The blue dotted line 
is an example of a flavor universal decay to two sleptons of 
different masses. The long-dashed green line is for two DM 
components decaying to muons and smuons. The experimen- 
tal data shown is similar to that in Figure 1. The two signal 
components and the background in the slepton case, as well 
as the background are shown as dotted lines. 



tons+sleptons via s — ► l^P and to sleptons via s — > 
^P. The lepton final states lead to a hard high energy 
feature while the slepton final state leads to a smoother 
one at lower energies. The purple short-dashed line of 
figure 2 shows the flavor universal decay of a scalar and a 
fermion of mass 1.5 TeV and equal abundance. Here the 
slepton masses are universal at 130 GeV and the LSP 
mass is 100 GeV. As shown in [1] a similar decay may 
yield a spectrum that is in remarkable agreement with 
the double feature observed by ATIC. Since the spec- 
trum of superpartners and the mass of dark matter set 
the scale of both spectral features (see equation (3)), a 
non-trivial cross check may be made once the LHC mea- 
sures the masses of sleptons and the LSP. 

Another possible way in which GUT induced dark 
matter decays can produce multiple spectral features is 
by various cascades of supersymmetric particles. For ex- 
ample, if dark matter decays to two sleptons of different 
mass (II and Ir) both of which decay to a neutralino 
LSP, the upper edges of the injection spectra in the two 
cascades may be sufficiently different (see equation (3)), 
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leading to multiple features. In figure 2 we show the 
spectrum produced by a flavor universal decay of dark 
matter into 145 GeV and and 230 GeV sleptons which 
subsequently decay to a 120 GeV LSP and a lepton (dot- 
ted blue line). The branching fraction into the heavier 
slepton is 20%. Again, since both spectral features are 
set by relations of the form (3) , this interpretation may 
be tested at the LHC. 

Alternatively, multiple features may be produced by 
two cascades of the same slepton via different neutralino 
states. For example, a similar spectrum (the dotted blue 
line in figure 2) may be produced by a 5 TeV scalar dark 
matter decaying to sleptons with a mass of 200 GeV. 
The slepton has two dominant decay channels into two 
different neutralino states with masses of 170 and 100 
GeV [42]. 

Spectral features can also arise from decays of two 
different dark matter particles with significantly different 
masses. Despite their different masses, the electron flux 
from the decays of these two particles can be compara- 
ble if their relic abundance is generated from the decays 
of another particle. For example, in the 50(10) model 
discussed earlier, the relic number density ng and n s of 
the singlets s and s are generated through the decays 
of the components (with mass m) of the 16 m that are 
charged under the standard model. With singlet masses 
rag and m s for s and s respectively, their relic number 

densities satisfy ^ ~ ( m ~ m ° ) . I n this model, the 

J n s \ m—m s J 1 

s and s can decay only when s develops a vev. Their 
respective dimension 6 decay rates Tg and T s scale as 

~ I™*] . The relative electron flux from the two 

decays is ~ (^) 3 (™) S '. This ratio is O (1) 

for ms + 0.4 ra s J$ m J$ mj + 2.2 m s when rag 3> m s . 
For TeV scale m s and mg, this results in comparable 
electron fluxes for ra within a TeV of mg + m s . The 
observation of these features in the electron spectrum 
can thus lead to measurement of SUSY parameters and 
GUT physics. In figure 2, we show the spectrum from 
the decays of a heavy 5 TeV scalar, s, and a 660 GeV 
fermion, s, into s — > p p^ and s — > A i± A =F - The decay 
rates are Tg = 8 x 10~ 27 s" 1 and T s = 3 x 10~ 27 s" 1 . 
The smuon was taken to be at 200 GeV and the LSP at 
90 GeV. In this case, the heavier components of the 16 m 
are around ^5.5 TeV. 

All of the scenarios described above are consistent 
with the qualitative shape of the PAMELA excess. The 
positron fraction for some of these cases is shown in fig- 
ure 3. Before proceeding to photon signals we will dis- 
cuss technical aspects of the figures above. We used 
GALPROP [6] for generating backgrounds and for prop- 
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FIG. 3: The positron fraction produced by some of scenar- 
ios discussed in the text - decays to smuons (Scenario 1, solid 
black), a universal lepton-slepton decay (Scenario 2, short- 
dashed purple) and decays to sleptons (Scenario 2, dotted 
black) . The background positron fraction from two propaga- 
tion models is shown in gray, thin dotted lines. The PAMELA 
data is shown (black circles). 



agation of the dark matter signal. In figure 1, we as- 
sumed the convective diffusion propagation model (DC) 
of [7]. In figure 2, a harder propagation model was used 
throughout (model B of [8]). The slope of the back- 
ground electron flux, which may significantly affect both 
the positron fraction and the total flux was chosen to be 
~ —3.3 at 20 GeV, in agreement with observations and 
the large uncertainties [9]. DarkSUSY [10] was used for 
producing injection spectra. 



Diffuse Gamma-rays 

Any decay scenario that produces charged parti- 
cles must produce gamma-rays from final state radia- 
tion (FSR) off those charged particles. Fig. 4 shows the 
gamma ray spectra from some of the models discussed 
above. The shape of the FSR spectrum is directly re- 
lated to the shape of the primary charged particle injec- 
tion spectrum from dark matter decay [1]. We do not 
include gamma-rays from inverse-compton scattering of 
starlight off the high energy electrons and positrons from 
the dark matter decay. This does not usually give as hard 
a spectrum as FSR, and presumably falls off faster than 
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FIG. 4: The gamma-ray spectra from final state radiation 
and r decay. The solid (black) and dot-dashed (red) lines 
are as in Fig. 1. The dashed (green) and dotted (blue) lines 
are as in Fig. 2. Only the signals are plotted, not signal 
plus background. Gray is the expected background from [12]. 
These are shown at a galactic latitude b — 60° and longitude 
I = 0°. 



FSR with galactic latitude as the density of starlight de- 
creases off the plane of the galaxy (but see [13]). This 
contribution is difficult to calculate off the plane of the 
galaxy due to the anisotropic flux of starlight [13, 14]. If 
a decay mode contains r's we do include the full spec- 
trum from these, e.g. the photons from 7T°'s produced in 
the r decay generated using [10]. 

This diffuse gamma-ray signal can determine 
whether the observed electron/positron excesses do in- 
deed arise from dark matter. Dark matter decays give 
rise to an isotropic gamma-ray signal with a hard, high- 
energy spectral feature such as an edge coming from 
FSR. The shape of the spectrum is the same everywhere 
across the sky. The intensity varies slightly, with a de- 
pendence on the angle of observation that is determined 
by the known density of dark matter in the galactic halo 
and so is only uncertain at the galactic center. If the 
electron/positron excesses arise from dark matter decay, 
the gamma-ray signal is probably strong enough to be 
observable at Fermi [1]. Observation of such a signal 
would be impossible to explain by any known astrophys- 
ical mechanism other than dark matter. 

These gamma-ray observations may even help dis- 
tinguish decays from annihilations, since decays produce 
a more isotropic flux than that from annihilations. Thus, 
the gamma-ray signal from decaying dark matter can be 
observed cleanly with measurements at high galactic lat- 
itude, off the plane of the galaxy, where there is little 
astrophysical background. 



The gamma-ray spectrum carries important infor- 
mation about the nature of the DM decay that is ob- 
scured by the electron/positron spectra. For example, 
the two scenarios in Fig. 1 produce almost identical elec- 
tron/positron spectra even though the underlying high 
energy physics is different. From Fig. 4 we see that these 
two are easily distinguished by their gamma-ray spectra. 
Further, the shape of the FSR spectrum can provide a 
measurement of the mass of the dark matter particle 
and potentially the masses of its decay products such 
as superpartners. For example, structure in the elec- 
tron/positron spectrum in Fig. 2 appears as well in the 
gamma-ray spectrum as in Fig. 4. In fact, the gamma- 
ray spectrum can provide extra information since FSR 
photons can come from any interior line in the decay 
chain that is on-shell. For example, in the dashed (green) 
curves in Figs. 2 and 4 a 300 GeV muon is produced 
which is directly visible as an edge in the gamma-ray 
spectrum. Thus, the gamma-ray spectrum can provide 
a probe of TeV scale physics complementary to the elec- 
tron/positron spectra. 

Neutrinos 

When the Dark Matter particle decays to charged 
leptons, neutrinos will also be produced due to either 
SU(2) invariance or subsequent decays of the produced 
muons or taus. In Fig. 5, we present the flux in Ice- 
Cube [15] for most of the above examples that predict 
neutrinos by SU(2) invariance. The solid black (dotted 
blue) line comes from s — > ^ M ( e ) + ^(e)i while the dot- 
dashed red line comes from s — > v r + v r . Because the 
neutrinos are produced at galactic distances, the flavor 
ratios on the earth are 1:1:1. We do not include the 
flux produced by charged lepton decays as they are sub- 
dominant. We take a bin size of 0.2 in Log 10 (of^) as 
it appears in [15], but the energy resolution could be 
as low as 0.4 in Log 10 (^) [16]. In the 100 GeV to 
several TeV range there is a large atmospheric neutrino 
background that drops rapidly with energy as a power 
law, E~ 3 . For example, taking into account the effective 
area of IceCube [16] and integrating over a bin of size 
0.4 in Log 10 (ofy) centered around the neutrino energy, 
5a discovery of a one TeV neutrino line could be possi- 
ble within roughly a year of observation time, when the 
Dark Matter lifetime is 10 26 sec. In this estimate we 
have not included systematics in the measurement pro- 
cess. Spectral information or distinguishing between dif- 
ferent scenarios is going to be harder to deduce because 
of the atmospheric neutrino background and IceCube's 
energy resolution. It is worth noting though that both 
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FIG. 5: The average z/ M flux produced from a Dark Matter 
particle decay in different examples. The solid (black) and 
dot-dashed (red) lines are as in Fig. 1. The dotted (blue) 
line is as in Fig. 2. AMANDA data is shown in gray [15]. 
The bin size is taken to be 0.2 in Log 10 (ofy), as it appears 
in [15]. 



the IceCube energy resolution as well as the atmospheric 
background subtraction could be greatly improved in the 
near future, increasing the potential of the experiment 
to observe and study astrophysical signals of TeV scale 
Dark Matter. 



CONCLUSIONS 

The decaying dark matter scenarios explored in this 
paper can naturally explain the cosmic ray spectra ob- 
served at HESS, PAMELA and ATIC without the need 
for astrophysical boost factors or additional new energy 
scales [1, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 
30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40]. In particular, 
the softer electron spectrum observed by HESS can fit 
naturally with the observations of PAMELA. Substruc- 
tures in the electron spectrum, correlated with substruc- 
tures in the photon spectrum, are also possible in these 
theories. Such substructure may also be correlated with 
the spectrum of superpartners observed at the LHC. Ob- 
servation of these substructures could open a window 
into the spectrum of TeV mass particles. The HESS ob- 
servation suggests a heavy (> TeV) DM mass, too heavy 
to be produced at the LHC. In this case, astrophysical 
observations would provide our only probe of dark mat- 
ter at such a high energy scale. 

The same theories typically also contain particles 
decaying during big bang nucleosynthesis through di- 
mension 5 operators with lifetime r ~ 8n T m 7 s. 
Such decays are recorded by a change in the primordial 



light clement abundances and may explain the anoma- 
lous observed Li abundances, opening another window 
to unification [1]. 
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It is generally interesting to point out that the HESS 
shape agrees well with dark matter decaying or annihi- 
lating to taus. 

It is unclear to what extent a hadronic decay of the next- 
to -lightest neutralino is in tension with the PAMELA 
anti- proton measurement [11]. The next-to-lightest neu- 
tralino may also decay to taus without producing an- 
tiprotons if the lightest stau is accesible. 



